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Abstract Analytic expressions for the speed, flux, microrotation, stress, and couple stress in a 
micropolar fluid exhibiting steady, symmetric and one-dimensional electro-osmotic flow in a uniform 
cylindrical microcapillary were derived under the constraint of the Debye-Hiickel approximation, 
which is applicable when the cross-sectional radius of the microcapillary exceeds the Debye length, 
provided that the zeta potential is sufficiently small in magnitude. As the aciculate particles in a 
micropolar fluid can rotate without translation, micropolarity influences fluid speed, fluid flux, and 
one of the two non-zero components of the stress tensor. The axial speed in a micropolar fluid 
intensifies as the radius increases. The stress tensor is confined to the region near the wall of the 
microcapillary but the couple stress tensor is uniform across the cross-section. 
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1 Introduction 

Electro-osmotic flows on the micrometer scale occur in many technoscicntiflc settings, including: 
microchanncls in chips for protein analysis and intravenous drug delivery [1]; microchanncls in bio- 
logical and chemical instruments [2]; micropumps, microturbincs, and micromachincs [3-5]; injectors 
of detoxification agents [6] ; and water desalinators [7] . These fiows occur by virtue of the movement 
of free ions in the liquid in the electric double layer (EDL) when an electric field is applied, and, in 
turn, cause bulk motion of the entire fiuid. The EDL is an infinitesimal region between a layer of 
charges of one polarity on the electrolytic-liquid side of the solid- liquid interface (i.e., the wall of the 
channel) and a layer of charges of the opposite polarity on the solid side. 

The discovery of electro-osmosis is just about two centuries old [8]. Fifty years after that event, 
electro-osmotic flow in a channel was experimentally found to be proportional to the applied current, 
which prompted Helmholtz to develop the EDL theory in 1879. Two decades later, the EDL thickness 
was measured as much smaller than the cross-sectional dimensions of the channel. In 1923, Debye and 
Hiickel published their analysis on the distribution of ions in a low-ionic-energy solution by using the 
Boltzmann distribution for the ionic energy [9], a seminal work that is widely used even nowadays. 

Many fluids are not simple Newtonian fluids — wherein stress is a symmetric tensor of the second 
rank. Instead, these fluids are micropolar as they contain aciculate particles that can rotate about 
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axes passing through their centroids [10,11]. Therefore, not only is stress asymmetric in micropolar 
fluids but they also sustain body couples. 

Micropolar fluids are exemplified by colloidal suspensions, liquid crystals, and epoxies [12]. Blood 
too is micropolar [13,14], along with other body fiuids containing particulate materials. As body 
fluids are subjected to electric fields in labs-on-a-chip [15], electro-osmotic fiows of micropolar fiuids 
in microchannels [16,17] are of emerging technoscientific importance in nanomedicine [18,19]. 

Our present interest lies in the spatial characteristics of fiuid speed, stress, microrotation, and 
couple stress in a micropolar fiuid flowing steadily in a microcapillary of circular cross-section with 
a cross-sectional radius exceeding the Debye length. With the assumption that zeta potential due to 
the EDL is sufficiently small, the Debye-Hiickel approximation can be applied to obtain analytical 
results. Analysis of these results would clearly show how the micropolarity infiuences fluid flow. 

This paper is organized as follows: The formulation of the relevant boundary- value problem is 
presented in Sec. 2, while Sec. 3 contains the description of analytical solution of the problem based 
on the Debye-Hiickel approximation. Section 4 contains numerical illustrations of the obtained 
analytical results and discussions thereon. The main conclusions are summarized in Sec. 5. 



2 Basic Analysis 

In a micropolar fluid, the stress tensor a' is deflned as a function of position r' as 

^'(r') = -p'{r')l + AIV • V'(r') + + x/2) { V'V'(r') + [V'V'(r')]^} 

-(x/2Ux [V'x V'(r')]+Xlxv'(r'), (1) 
and the couple stress tensor m' as 

^'(r') = alV ' v'(r') + /?V'v'(r') + 7 [V'v'(r')]^ . (2) 

Here, / = V'r' is the idempotent; p' is the hydrostatic pressure; V and v', respectively, are the fluid 
velocity and the microrotation; a, /3, and 7 are the three spin-gradient viscosity coefficients; and fi 
and X, respectively, are the Newtonian shear viscosity coefficient and the vortex viscosity coefficient 
related by the inequality 2/i -I- x > 0, where x ^ [12, p. 14]. The parameters v', x, Oi, P, and 7 are 
null- valued in a simple Newtonian fluid. The superscript T denotes the transpose. 

We examine here the steady flow of a micropolar fluid in an inflnitely extended microcapillary 
of cross-sectional radius R, such that its axis coincides with the z' axis of the cylindrical polar 
coordinate system r' = {p',0',z'). The study is undertaken under the following assumptions: (i) 
the zeta potential is uniform in the microcapillary; (ii) the surface of the microcapillary is perfectly 
insulated and impermeable; (iii) the applied electric field is spatiotemporally uniform electric field 
and aligned parallel to the axis of the microcapillary; (iv) the fluid is ionized, incompressible, and 
viscous; (v) the flow is fully developed, steady, laminar, axial, and radially symmetric; (vi) the effect 
of gravity is negligible; (vii) neither a body couple nor a pressure gradient is present; (viii) the Joule 
heating effects are small enough to be ignored; and (ix) R is much greater than the Debye length A d ■ 

Under these conditions, our starting point comprises the following three equations of micropolar- 
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fluid flow [12]: 

V'-V'(r')-0, (3) 

V • g'(r') + pl(r') Kpp = P>n[V'(r') • V']V'(r') , (4) 

V • ^{r')+L X ^'(r') = p„jo[V'(r') • V']v'(r') . (5) 
Introducing Eqs. (1) and (2) in these three equations, we get 

V'-V'(r')=0, (6) 



+ x) V X [V X V'(r')] + xV X v'(r') + p^(r') E^^^ 

= p„,[V'(r') • V']V'(r'), (7) 



(a + /3 + 7)V' [V • v'(r')] - 7V' x [V x v'(r')] - 2xv'(r') 

+XV' X V'(r') = p„jo[V'(r') • V']v'(r') . (8) 

Here, Ej^^p is the apphed electric field, whereas pm and jo, respectively, are the mass density and the 
microinertia. The microincrtia is null-valued in a simple Newtonian fluid. 

In the absence of a significant convective or electrophoretic disturbance to the EDL, the charge 
density /Ce(r') is described by a Boltzmann distribution, and takes the following form for a symmetric, 
dilute, and univalent electrolyte [20]: 

p'^{r') = ~2zoeno sinh [zoei^'{r')/kBT] . (9) 

Here, Zg is the absolute value of the ionic valence, ip' is the electric potential, e is the charge of an 
electron, Ug is the number density of ions in the fluid far away from any charged surface, fc^ is the 
Boltzmann constant, and T is the temperature. With e denoting the static permittivity of the fluid, 
the charge density and the electric potential are also related by the Gauss law 

V'V(r') = -Pe(r')A; (10) 

thus, 

V'V(r') = {2zoeno/e) sinh [zoeib'{r')/kBT] . (11) 

The Debye length \o = {zoe)~^ (eA:BT/2no)^^^ is assumed in this paper to be much smaller than 
the cross-sectional radius R. Furthermore, we set Ej,pp = z' Eq as the applied electric field. 

Before proceeding, let us define the non-dimensionalized quantities 

r = r'/i?, Y = V'/U, w = {R/U)V ^ 

^= (i?/;7)(/i + x)"y , m={R^/lU)^ \. (12) 

^ = V''/V'o, pe = (i?Ve^o)pL 

Here, r = {p, 9, z) with p = p' / R, 9 = 9', and z = z' /R; 

U= — 13 

Ai + X 
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is a characteristic speed [16]; and ipo is called the zeta potential [20] which is assumed to be temporally 
constant and spatially uniform at the wall p = 1 of the microcapillary. With these quantities, Eqs.(6)- 
(8) and (11), respectively, simplify to 



V • V(r) 0, 

-V X [V X V(r)] + fciV X v(r) - (r) z' 

= i?e[V(r) . V]V(r), 

-V X [V X v(r)] - 2/s2v(r) + fcgV [V • v(r)] 
= i?o[V(r) • V]v(r), 

VV (r) = {ml/ao) sinh [ooV (r)] , 



fczV X V(r) 



where 



fci 

Re 



X 

M+X ' 



k2 

Rn 



xEL 

7 ' 

1 



m„ 
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(14) 

(15) 

(16) 
(17) 

(18) 



Here, fci couples the two viscosity coefBcients, ^2 and are normalized micropolar parameters. Re 
may be called the Reynolds number, Ro may be called the microrotation Reynolds number [12], and 
ao is the ionic-energy parameter [20]. The Gauss law (10) can now be written as 



V2^(r) = -p,(r). 



(19) 



The general condition of radial symmetry implies that d/d6 = 0. As the length of the microcapillary 
is infinite while its cross-sectional radius is finite, we set d/dz = 0; thus, V(r) = 'V{p) and v(r) = 
v(p) Since the flow is assumed to be axial and laminar as well, it follows that 'V{p) ~ Vj(/9)z'. 
Equation (14) is then automatically satisfied, whereas Eq. (17) reduces to 



d_ 
dp 



dp 



and Eq. (19) to 



Using Eq. (21), we obtain 



Pe{p) 



pml 



dp 



sinh [ttoV' (p)] 



d 
dp 



dp 



and 



dVM ^ d^jp) 

dp 

dvzjp) _ 
dp 



dp 



kive{p) 



from Eq. (15). After using Eq. (23), we get the following three equations from Eq. (16): 



d_ 
dp \ 

dp\ 

■"zip) 



p 



dp 



dp 



[pVp{p)]^~2k2Vp{p)^G., 

[pve{p)] \ - 2k2Ve{p) - k2 



dVzjp) 
dp 



(20) 
(21) 

(22) 
(23) 

(24) 

(25) 
(26) 
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According to Eq. (24), Vp is not influenced by the axial flow. Therefore, we set 



vM=^^ (27) 

and focus our attention on Eqs. (20), (22) and (25) involving ^{p), Vz{p) and vg{p). 

We need to provide appropriate conditions for these three quantities. First, by virtue of the 
definition of the zeta potential, 

V(l) = l; (28) 

next, the condition 

V'(O) < 1 (29) 
is engendered by the assumption R:^ Xd [8, p. 187]; finally, fiow symmetry dictates that 



#(P) 



= . (30) 

p=0 



dp 

The no-slip boundary condition [20, pp. 100-101] on the surface p = 1 of the capillary is 

K(l) = (31) 

and the symmetric flow requires that 

0. (32) 



dV^ip) 



dp 



p=0 



Finally, the boundary conditions imposed on microrotation are: 

veil) =l3o^^^ , veiO) = (33) 
dp p=i 

where /3o G [^IjO] is some constant. Although some researchers have ignored microrotation effects 
near a solid wall by setting (3o = [12,21], others have held that that (3o < because the existence of 
the boundary layer requires that the shear and couple stresses on a wall must be high in magnitude 
in comparison to locations elsewhere [22]. Furthermore, the limiting case of /3o = —1 accounts for 
turbulence near the wall [23]. As we think that /3o < may be reasonable if electro-osmosis occurs 
because of the presence of the EDL near the wall p = 1, results are provided in this paper for /3o 7^ 0. 



3 Debye— Hiickel Approximation 

If the electric potential energy is small as compared to the thermal energy of the ions, i.e., \zoeipo'4'\ ^ 
jfesTj, then [ckoV'l ^ 1; accordingly, 

sinh{aoip) ~ ctoip j (34) 



which is called the Debye-Hiickel approximation [20] . In light of this approximation Eq. (20) simplifies 
to 

P , 2 + ~i — -rn^P^iP), (35) 



whose solution 

loiruop) , . 

^KP) ~ -Ti — r ■ (36) 
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satisfies the boundary eonditions (28) and (30); here and hereafter, /«( • ) is the modified Bessel 
function of the first kind and order n [24]. As nio ~ R/Xd oo, we see that "0(0) « l/Ioinio) — > 0, 
thereby fulfilling the requirement (29). 

Integrating both sides of Eq. (22) with respect to p and using the boundary conditions (30), (32) 
and (33)2, we obtain 

— = -kive{p) — . (37) 

dp dp 

Eliminating dVz/dp from Eq. (25), we next get 



^.dMp) ^ _ ^ (2 - kM]Mp) + fc.P^^ - . (38) 

ap^ dp dp 



Finally, using the Debye-Hiickel solution (36) for ip{p), we reduce Eq. (38) to 

P , 2 + P—i (1 + Kp ) ve{p) = — p7 ^fc2"ioP , (39 

dp^ dp loymo) 

where 

fco = v/(2 - fci)fc2 . (40) 

The solution of the homogeneous counterpart of Eq. (39) is ci/i(fcop) + C2ii'i(fcop), where Kn{ ' ) 
is the modified Bessel function of the second kind and order n [24] , while Ci and C2 are constants to 
be determined later. The method of variation of parameters then yields the particular solution of 
Eq. (39) as k2mo (fcg - ml)'^ /i(mop)//o(mo), after the identity J„(^)i<'„+i(^)+/„+i(^)X„(^) = 
[1, Eq. 9.6.15] has been exploited. Since Ki{(,) — >■ 00 as ^ — >■ cxd but /i(0) = 0, satisfaction of the 
boundary condition (33)2 requires that C2 = 0. Therefore, the complete solution of Eq. (39) is 



ve{p) = cihikop) + -2 2 -Tf T ' ^^^) 



Substitution of Eqs. (37) and (41) in the boundary condition (33)i leads to 

Finally, using Eqs. (36) and (41) in (37), and exploiting the boundary condition (31), we get 
V.{p) = ci-i [loiko) ~ loikop)] + 1 + 

Kq \ Kq it 

which automatically satisfies the requirement (32). 



(43) 



Thus. Eqs. (41)-(43) constitute the solution of the boundary-value problem for steady fiow of a 
mieropolar fiuid when the Debye-Hiickel approximation (34) holds, for all values of /3o. 

On setting x = for a simple Newtonian fiuid, we get fci = fc2 = = 0. Accordingly, Eqs. (41) 
and Eq. (43) simplify to 

.«(p)^0, K(p) = l™:|^, (44) 
which is a known result [20, p. 102]. 
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4 RESULTS AND DISCUSSION 



All calculations were made with the following material properties fixed: Uq 

= 1, e = lOeo, Co = 8.854 x fQ-^^ p ^-i^ /i = 3 x 10"^ and 7 = fO"'' kg m s"^ The 

temperature was fixed at T = 290 K. Since the Boltzmann constant ks = 1-38 x fO^^'^ J K^^ 
and the electron charge e ~ 1.6 x 10~^^ C, the Dcbye length Xd = 10.72 nm. Consistently with 
the assumption that R ^ Xd, wc chose R £ [107.2,5360] nm so that nio G [10,500]. We fixed 
ipo ~ —25 X 10^'^ V, which is about the upper limit for the Debye-Hiickel approximation to be valid 
at about the room temperature [25, p. 25]. The parameters (3 € [—7,7] and ki G [0,0.95] were kept 
as variables, after noting that /si — > 1 as x — >■ 00. The magnitude of the applied electric field was 
fixed at Eo = 10^ V m^^, which is a reasonable practical value. The dependences of the relevant 
components of the fluid velocity, microrotation, stress tensor, and couple stress tensor on TOq, ki, 
and Po were investigated for steady flow. 

Besides the speed V^{p') and the microrotation Vg{p'), the only non-zero components of the stress 
tensor and the couple stress tensor, according to Eqs. (1) and (2) are 



The non-dimensionalized form of these equations are 

a,p(p) = (l-fci)^-Me(p) 



dp 



(45) 



(46) 



where kj = 13/ j G [—1, 1] 



4.1 Fluid speed 

In order to examine the infiuence of parameters ttIq, fci, and /3o on the fluid speed, we computed 
V^{p) using its Debye-Hiickcl expression (43). Some representative plots of V^{p) vs. p are provided 
in Figure 1 for e {50, 500}, fci e {0, 0.5, 0.95}, and € {-0.01, -0.1, -0.5, -1}. 

Figure 1 indicates that the speed of a simple Newtonian fluid (fci = 0) depends on nio but not on 
/3o, in conformity with Eq. (44)2. Moreover, since < 7^7^ < 1 for all e [50, 500] and p G [0, 1], 
< V^{p) < U hy virtue of the same equation. In addition, for any p g [0, 1), V^{p) increases with 
nio and the maximum fluid speed exists in the center (p = 0) of the microcapillary. 

Figure 1 also shows that the speed of a micropolar fluid (fci 0) decreases as ki increases for 
all Too but not for all /3o. In addition to the fluid speed increasing with TOq (just as for a simple 
Newtonian fluid), the fluid speed also increases with \(3o\ for all nio ^ 1 and fci > 0. 
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Figure 1: Variation of V^{p) with p when (a-d) ttiq = 50 or (e-h) irio = 500, for fci = (solid curves), 
ki = 0.5 (dotted curves), and fci = 0.95 (dashed curves), (a, e) /3o = —0.01, (b, f) (3o = —0.1, (c, g) 
/3o = -0.95, (d, h) = -1. 



The foregoing trends are in accord with the results of an analytical investigation in the central 
portion of the microcapillary. With the assumption that mo > 10, Eq. (43) yields 

VM^u(l + kJ-^-^-^(-^ + -^) I'-^f^]]) ■ (47) 
V {K~ml ko \kl-ml 1 + fci^o/ L /i(fco) \]J 

For simple Newtonian fluids (fci = 0), Eq. (47) shows that V!^{0) = U is independent of ruo as well 
as of, as expected, jSo- Plots of V!.{Q) vs. rUo in Figure 2 indicate that the fluid speed at the center 
of the microcapillary increases with TOq, and therefore with i?, for micropolar fluids (fci ^ 0). When 
microrotation effects can be neglected at the wall of the microcapillary — i.e., when = — Eq. (47) 
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simplifies to 



K'(o) 



/3o=0 



U I 



loih) - 1 



/i(fco) 



(48) 



and Figure 2 shows that the effect of TOq on V^{Q) is very weak but the effect of fci is considerable. 
As |/3o| increases, ^^'(0) acquires a stronger tendency in a micropolar fluid to increase with mo, which 
can be concluded from Figure 2. 



(a) 



(b) 



(c) 



Figure 2: Variation of V^{0) with TOq for fci = (sohd curves), fci = 0.5 (dotted curves), and fci 
(dashed curves), (a) /3o = 0, (b) f3o = -0.2, (c) = -1- 



0.95 



The fluid-speed gradient dV^/dp is not only an important factor in microrotation per Eqs. (33) i 
and (37), but it also influences several components of the stress tensor identified in Eqs. (46). The 
Debyc-Hiickcl expression for dV^/dp follows from Eqs. (36), (37), and (41) as 



dp 



cikiIi{kop) + ?7io 1 + 



(49) 



Now, this gradient vanishes as p — >■ 0, i.e., at the centre of the microcapillary, in conformity with 
boundary condition (32). But it has high magnitudes near the wall, as can be gathered by setting 
p = 1 in Eq. (49) to obtain 



dV^{p) 



dp 



Uvrio h{mo) 



1 + kiPo h{mo) 



When mo > 10, the foregoing expression can be further approximated as 



dV^ip) 



dp 



p=i 



Unio 
1 + kiPo 



(50) 



(51) 



evincing a direct proportionality with nio (or R) for all |/3o| — for both simple Newtonian and microp- 
olar fluids. 



4.2 Fluid flux 

From an engineering perspective, the fluid flux should be considered in addition to fluid speed. In 
the present context, it is defined as 

Q' = 2^/ p'V;{p')dp' ^2ttR^U [ pV,{p)dp. (52) 
Jo Jo 
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On substituting the Dcbye-Hiickel expression (43) in the integrand on the right side of Eq. (52), we 
get 



which simplifies to 



Q' ~ ttR^U 



ko 

Clfcl 



1 



kik2 



hiko) 



kik2 



h{mo) 

2 

nic 



(53) 



(54) 



for nio > 10. Clearly then, the fluid flux depends on fci, /3o, and mo- 

The dependence of Q' on Po can be identified using Eqs. (42) and (53). We can write Q' as the 
sum of two parts, one of which is independent of Po and the other depends linearly on Po/ (1 + fci/3o)- 




400 500 



SOOOi ^ 

7000 



^ 5000 
E 



; 3000 

r 2000 



200 300 400 500 



Figure 3: Variation of Q' with TOq for fci 
(dashed curves), (a) /3o = —0.01, (b) /3o = 



(solid curves), fci 
-0.2, (c) Po = -1. 



0.5 (dotted curves), and fci — 0.95 



Although for simple Newtonian fluids (fci = 0), Eq. (54) yields 



1 



Q' 



ttR^U 



(55) 



for nio 3> 1, the relationship of Q' and nio is more complicated when fci ^ 0. The variation of Q' 
with respect to mo is illustrated in Figure 3 for nine different combinations of fci and Po- This figure 
shows that Q' increases with nio (and, therefore, with R) for simple Newtonian fluids; the same trend 
exists for micropolar fluids, regardless of the value of Po G [—1,0]. Furthermore, Q' intensifies with 
increasing \Po\ for all nio 6 [50, 500] and fci G (0, 1), the intensification rate dQ' /d\Po\ itself increasing 
concurrently. 

Figure 4 shows the dependency of Q' on fci for 12 different combinations of Po and rrio. Above 
a threshold value of \Po\ which is quite small, the following trend is followed: from the value Q'q ~ 
nR^U{l — 2/mo) at fci = 0, the fluid flux Q' increases linearly with fci for fci <C 1, then increases 
nonlinearly with fci to a maximum value Q' at fci = fci, and finally drops monotonically to Q'l at 
fci = 1, for Po G (—1,0] and nio 3> 1. Whereas Q' increases with both nio and \Po\, fci is almost 
totally independent of nio but increases with \Po\. Calculations show that Q[^ni2, with the exponent 
n lying between 2 and 3 and increasing with \po\- Furthermore, for Po = —1, Q' increases linearly 
with fci for fci G [0, 1] and fci = 1. 
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Figure 4: Variation of Q' with fci for Po = —0.01 (solid curves), (3o ~ —0.2 (dotted curves). Po = —0.99 
(dashed-dotted curves), and Po ~ —1.0 (dashed curves). (a) irio = 50, (b) TOq = 275, and (c) TOq — 500. 

4.3 Microrotation 

The microrotation v'g{p) exists in micropolar fluids but not in simple Newtonian fluids. Accordingly, 
representative plots of v'g vs. p, computed using the Debye-Hiickel expression (41), are presented in 
Figure 5 for TOq 3> 1 and ki > 0. This figure indicates that v'g{p) 

(i) increases linearly with p, 

(ii) increases with |/3o| for all fci £ (0, 1) and p G (0, 1], and 

(iii) decreases as ki increases for aU nio £ [50,500] and /3o G [—1,0). 

By virtue of the boundary condition (88)2, the microrotation is weak in the the central part of 
the microcapillary, but it is maximum at the wall p = 1 for all ttIq 3> 1, /3o £ [~liO)i s-^d fci > 0. 
Indeed, Eq. (41) yields 

p=l XD{l + kiPo) 

This equation shows that the microrotation at the wall of the microcapillary is independent of i?f, 
which is in agreement with Figure 5. This figure also shows that, for /3o = — 1, at p = 1 is 
independent of fci, which agrees with the result 

v'eif) ^ (57) 

obtained after using the definitions (13) and (18) in Eq. (56). 

All of the foregoing conclusions about the microrotation can be encapsulated as the approximation 

it 1 + fciPo 

which is valid for TOq 3> 1. Its derivation proceeds as follows: A comparison of magnitudes for large 
•nio suggests that fc2/ (fc'o — jtIq) can be neglected in favor of f5o/ (1 + fci/3o) while Ii{mo) / Io{mo) — 1 
in Eq. (42), leading to 

l + kiPo h{ko) 
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p 



p 



Figure 5: Variation of (f>v'g{p), where the normahzation factor (j) = (68|/3o|) ^, with p when (a-d) 
Too = 50 or (e-h) TOq = 500, for fci = 0.01 (sohd curves), fci = 0.5 (dotted curves), and ki — 0.95 
(dashed curves), (a, e) = -0.001, (b, f) /?„ = -0.1, (c, g) /3o = -0.99, (d, h) Po = -1- 



A similar argument leads to the neglect of the second term on the right side of Eq. (41) yielding 

l + kiPo hiko) 

Even for very large TOq, ko is very small so that /i(fco) ~ fco/2 and /i(fcop) ~ kQp/2; Eq. (58) then 
follows for Too ^ 1. 
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4.4 Stress tensor 



After the use of Eqs. (41) and (43) in Eqs. (45)i and (45)2, we obtain the non-zero components of 
the stress tensor as 



Xr 



loim-o) 



and 



eipoEo ( cifci (2 - ki) 



A 



D 



hikop) 



2 1 



(61) 



(62) 



consistently with the Debye-Hiickel approximation. For a simple Newtonian fluid (fci = 0), from the 
foregoing expressions we get (j'p^ip) = <^'zp{p)^ which conforms to the symmetry of the stress tensor 
in the absence of micropolarity. In a micropolar fluid (fci ^ 0), the stress tensor of Eq. (1) has both 
symmetric and skew-symmetric components, which explains why (j'pz{p) 7^ ^'zp{p)- 

What is really surprising is that cr'^^^p) is the same for simple Newtonian as well as micropolar 
fluids. This component of the stress tensor depends on m^, but neither on ki nor on Pq. Furthermore, 
per Eq. (61), it is absent on the axis of the microcapillary because /i(0) = 0. Indeed, as can be 
deduced from Fig. 6, cr'pzip) is negligible in most of the microcapillary and becomes significant in 
magnitude very near to and on the wall, with 



dp) 



An 



(63) 



As nio increases, the region in which the magnitude of a' {p) is significant shrinks. 




Figure 6: Variation of cjpzip) with p when (a) nio — 50 or (b) TOq ~ 500. There is no dependence on 
either fci or /So- 

The difference (j'pzip) — <^'zp{p) entirely due to micropolaity, as becomes evident on comparing 
the spatial profiles of cr'zpip) in Fig. 7 with the spatial profiles of cr'pzip) in Fig. 6. Clearly, a'^p{p) 
depends on ruo, ki, and /3o- Now, cr^p(O) = but 

' f \ _ ej^oEo (1 - fci - kiPo) ( 1 \ 

for nio :» 1, and Fig. 7 shows that a'^p{0) > cr'zpip) > <^'zpi^) fo'' P (0' !)■ The region around 
the central axis of the microcapillary can be called a cr^p-free zone, which shrinks rapidly as |/3o| 
increases. 

Figure 7 also shows that \(7'^p\ not only increases with ki for TOq 3> 1 but it also increases with 
nio ^ 1 for all fci G [0, 1] and /3o G (—1, 0). As Po ^ 0, the distinction between the micropolar and 
the simple Newtonian fluids disappears for a'^p(p). 
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(a) 



'■■ (b) 



(c) 



(d) 



(e) 




(h) 



Figure 7: Variation of cr'^p{p) with p when (a-d) nio = 50 or (e-h) rrio = 500, for fci = (solid curves), 
ki = 0.5 (dotted curves), and ki = 0.95 (dashed curves), (a, e) Po — —0.001, (b, f) /?o = —0.1, (c, g) 
Po = -0.5, (d, h) = -1. 



4.5 Couple stress tensor 

The couple stress arises in conjunction with the microrotation, as delineated by Eq. (2), and therefore 
cannot exist in a simple Newtonian fluid. Consistently with the approximations made for a micro- 
capillary, the only two non-zero components of the couple stress tensor are given by Eqs. (45)3,4 or. 
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equivalcntly, Eqs. (46)3.4. Substitution of the Dcbye-Hiickel approximation (41) for vg therein yields 



and 



+ fc8 



7C/ 



Cl 



fc? + 1 

kokjloikop) hikap) 



mok^Ioimop) ^ Ii{mop) 



(65) 



Cl 



k^ + 1 

kohikop) h{kop) 



mohiri-Lop) - ^^-^^himop) 
P 



where 



k^ - ml loinio) 



(66) 
(67) 



Equations (65) and (66) yield the difference 



7C/ 



"^pe(/5) - m'gp{p) = (^7 - 1) [cikoIo{kop) + ksmolo{mop)] 

between the two non-zero components of m'. This difference vanishes for all p G [0, 1] when /3 = 7 
(i.e., A;7 = 1), in agreement with Eq. (2). Furthermore, as 



m' (0) = -m' (0) 



jU_ 



(fcy - 1) (cifco + k^mo) 



(69) 



the condition /3 = 7 also implies that the couple stress tensor is nonexistent on the axis of the 
microcapillary. 

A more general but slightly approximate conclusion can be drawn from Eq. (58) for nio ^ 1 
as follows. That equation yields vg{p) ~ — /3o (1 + fci/3o)~ 'moP-, whose use in Eqs. (46)3,4 and (12) 
provides 



(70) 



for Too 3> 1. In other words, the couple stress tensor is skew-symmetric as well as uniform throughout 
the cross-section of the microcapillary. The approximate equations (70) were verified by directly using 
Eqs. (65) and (66) to plot the variations of m' Jp) and rrig (p) with p in Figs. 8 and 9, respectively. 



5 CONCLUDING REMARKS 

We formulated the boundary-value problem of steady electro-osmotic flow of a micropolar fluid in 
a microcapillary whose length is much greater than its cross-sectional radius. Analytical solution of 
this boundary- value problem was obtained under the assumption that both the Debye length and the 
zeta potential arc sufficiently small in magntitudc so that the Debye-Hiickcl approximation could be 
used. 
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Figure 8: Variation of m^g(p) with p when fcy = 0.5 and (a,b) nio = 50 or (c,d) rUo = 500, for 
ki = 0.1 (solid curves), ki = 0.5 (dotted curves), and ki = 0.99 (dashed curves), (a, c) j3o = —0.01, 
(b, d) Po = -0.99. 
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Figure 9: Same as Fig. 8, except that m'g [p) is plotted against p. 



As the aciculate particles in a micropolar fluid can rotate without translation, micropolarity 
must influence fluid speed. The derived expressions and the calculated data indicate that the fluid 
speed depends significantly on the viscosity coupling parameter ki , which mediates the (micropolar) 
viscosity coefficient and Newtonian shear viscosity coefficient. The axial speed of a micropolar fluid is 
below the speed of a simple Newtonian fluid for small values of the micropolar boundary parameter 
Po — which relates the velocity gradient and microrotation at the wall of the microcapillary — but 
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exceeds the latter for j3o close to — 1 . In addition, the axial speed is independent of the radius of the 
microcapillary when the fluid if of simple Newtonian type but not when it is micropolar, provided 
the Debyc length is fixed; indeed, the magnitude of the axial speed in a micropolar fluid intensifies 
as the radius increases. 

The fluid flux increases as either the microcapillary radius increases and/or the Debye length 
decreases, whether the fluid is simple Newtonian or micropolar — in the latter case, regardless of the 
value of the micropolar boundary parameter (3o. The flux of a micropolar fluid increases as the 
magnitude of that boundary parameter intensifies. 

Although microrotation greatly influences the speed and the flux, it vanishes on the axis of the 
microcapillary and increases linearly in the radial direction. Microtation increases with the magnitude 
of the micropolar boundary parameter j3o, but it decreases as the viscosity coupling parameter ki 
increases. Quite surprisingly, microrotation at the wall is independent of the cross-sectional radius 
of the microcapillary. Moreover, when the boundary layer is turbulent (i.e., f3o = —1), microrotation 
at the wall is also independent of ki. 

The stress tensor in the fluid has just two non-zero components, one of which is totally unaffected 
by the micropolarity of the fluid. That component does not exist on the axis and is largely confined to 
the region close to the wall. The other component is also absent on the axis and it gets progressively 
concentrated on the region close to the wall as (So 

Unlike all foregoing physical parameters, both non-zero components of the skew-symmetric couple 
stress tensor are uniform in a micropolar fiuid throughout the cross-section of the microcapillary. The 
couple stress tensor does not exist in a simple Newtonian fluid. 

Our conclusions are significant for the design of microcapillaries are the selection of materials 
for labs-on-a-chip. For instance, turbulence caused by mixing of a (micropolar) body fluid with a 
(simple Newtonian) reagent fluid is likely to result in higher electro-osmotically induced flux and 
axial speed in a microcapillary than if both fluids arc of the simple Newtonian type, suggesting 
that the microcapillary be designed with a larger cross-sectional diameter. Higher stress at the wall 
of a microcapillary transporting a micropolar fluid suggests that stiffer materials be used for the 
construction of labs-in-a-chip than if all fluids were to be simple Netwonian. As all of our conclusions 
apply when the zeta potential is sufficiently small in magnitude and the the cross-sectional radius of 
the microcapillary exceeds the Debye length, numerical solution of Eqs. (6)-(8) is required for more 
general situations. We plan to take up that investigation next. 
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